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System components and tunable parameters 



Program 


Function performed 




Parse BLAST output 


mark substitution heteropJasmy.pl 


Mark ambiguous base calls 




Extract a consensus sequence from a 
phrap polyphred run, replacing tow quality 
bases with *?\ and polyphred rank i 
bases with 'N r . 


computejcoverage.pl 


Calculates overall, forward strand, and 
reverse strand coverage for ail bases in 
the consensus sequence 




Count deletions in HV1 relative to the 
rCRS 


borderjndexpl 


Compute start and end positions of HV1 
and HV2 regions. 


fjx_mrtGtype_reponing_range-pi 


Outnuf a mtDNA profile based on input 
variants list and reportable range.. 


find bad traces from blast reportpl 


Reject reads that do not align 
appropriately to the rCRS 


determineRea3Types.pl 


Adda template name, template type, and 
primer type to phred output files. 


seq2delta_vs 


Align mtDNA profile to rCRS and report 
valiants according to nomenclature. 


caloilate_coverage_mitotype.pl 


Mask mtDNA consensus sequence based 
on required coverage and calculated 
coverage 


countjwl jnserts.pl 


Count insertions in HV1 relative to the 


generate_hv_mask_fastajiles.pl 


Extract the HV1 portion and HV2 portion 
of the mtDNA consensus sequence, 
based on the computed HV1 and HV2 
regions. 


flip_fasta.pl 


Invert an X masked sequence output by 
cross match 
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Assessment I Assessment Z 



Total Tme Positives: N at ^™^= Naunmj 


.16 


10 


Total True Negatives: A&to&kA- Acuaiiai 


77358 


: 110,354 


Total False Positives: N«mwwted~ Aaam*! 


95 


807 


Total False Negatives: A**^^ ~ Noum*! 


6* 


14 


Total Inconect 


0 


0 


Sensitivity: TP/(TP+FN) 


7173% 


41.67% 


Specificity: TW(TN+FP) 


99.88% 


99-27% 


Positive Predictive Value: TP/(TP+FP) 


14.41% 


. 1.22% 


Negative Predictive Value: TN/(TN+FN) 


99.99% 


99.99% 



1 Predominate base correctly cal led 

Differences take into account the consensus of the two analysts 
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